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Abstract 

The <?-function was introduced by Affleck and Ludwig as a measure of the ground state 
degeneracy of a conformal boundary condition. We consider this function for pertur- 
bations of the conformal Yang-Lee model by bulk and boundary fields using conformal 
perturbation theory (CPT), the truncated conformal space approach (TCSA) and the 
thermodynamic Bethe Ansatz (TBA). We find that the TBA equations derived by LeClair 
et al describe the massless boundary flows, up to an overall constant, but are incorrect 
when one considers a simultaneous bulk perturbation; however the TBA equations do cor- 
rectly give the 'non-universal' linear term in the massive case, and the ratio of (/-functions 
for different boundary conditions is also correctly produced. This ratio is related to the 
Y-system of the Yang-Lee model and by comparing the perturbative expansions of the 
Y-system and of the g-functions we obtain the exact relation between the UV and IR 
parameters of the massless perturbed boundary model. 
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1 Introduction 



Consider the partition function of a classical statistical-mechanical system defined on a 
cylinder of length R and circumference L. Among the characteristics of the model might 
be a bulk mass scale M and boundary scales depending on the boundary conditions a and 
(3 imposed at the two ends of the cylinder; we will highlight the role of these quantities by 
denoting the partition function Z a @(M, R, L). If R is taken to infinity with all other variables 
held fixed, then 

Z a0 (M,R,L) ~ A Qf3 (M,L)e- RE o^ M ^ , (1.1) 

where Eq 1tc (M,L) is the ground state energy of the model on a circle of circumference L. 
To derive this asymptotic i?-dependence, it is sufficient to treat the boundary conditions as 
boundary states | a ) ||] in a formalism where time runs along the length of the cylinder, and 
states are propagated by a bulk Hamiltonian H c \ rc (M,L): 

Z aP (M, R,L) = (a \ exp(-RH ciTC (M, L))\(3) . (1.2) 

At large R the contribution of the ground state | £1 ) dominates, establishing ( |1.1| ) and also 
giving 



The inner products appearing in ( |1.3| ) should in general contain a term corresponding to a 
free-energy per unit length, i.e. 

log( ( ^"i/2 ) = ~ Lfa + lo g(^(M, L)) . (1.4) 

This linear term can in principle be extracted unambiguously from the large L behaviour of 
log( ( O | a ) ). The question is then whether the functions log(<? a (M, L)) now contain universal 
information. 

In the case M = 0, i.e. for critical bulk, Affleck and Ludwig M pointed out that the UV 
and IR limits of these functions, log(g Q (0, 0)) and log(g Q (0, oo)), play the role of a generalised 
ground state degeneracy for the UV and IR conformal boundary conditions respectively. 
These can be easily calculated for many conformal field theories, and can enable one to 
identify the boundary conditions uniquely. 

The universality of the ^-functions is however a somewhat delicate issue when the model 
has a mass scale, either in the bulk or at the boundary. For example, one could imagine that 
in some calculational schemes the boundaries acquire a finite thickness, and so the effective 
cylinder length would decrease by some finite amount 5 to R — 5, in which case log(g Q ) would 
be altered 

log( 9a (M,L)) -> log(g a (M,L)) + ^c(ML) , (1.5) 
where c(ML) is related to the ground state energy on the circle by 

E™ C (M, L) = f hulk L - ^c(ML) . (1.6) 
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More generally, the fact that the ^-functions have to be extracted from subleading contri- 
butions to log Z makes their calculation much trickier than that of the corresponding bulk 
quantity c(ML). 

If the model is integrable, one might hope to be able to address these issues in greater 
depth. This was one of the topics of some work by LeClair et al. ||. Using the thermodynamic 
Bethe ansatz (TBA) technique, they proposed equations for log(<?), but up to now these have 
not been checked in detail. The main purpose of this paper is to provide some such checks 
in a concrete example, namely the scaling Yang-Lee model, making use both of conformal 
perturbation theory and of a recently-introduced generalisation of the TCSA to boundary 
situations Q. We find that the TBA ^-functions describes the massless flow well (up to 
an overall constant) but the massive flows are not well described. However, we find that the 
TBA (/-functions of |J do appear to describe correctly the ratio of the ^-functions for different 
boundary conditions of the same bulk model. Further, the TBA expression for this ratio has 
a very simple form in terms of the Y-function of the Yang-Lee model. This enables us to (a) 
provide a perturbative expansion for the massive Y-function and (b) find the coefficient in 
the relation between UV and IR parameters in the boundary Yang-Lee model model exactly, 
a result announced in 

We first discuss the conformal perturbation theory (CPT) expressions for log g and then 
how the truncated conformal space approach (TCSA) can be used to investigate these func- 
tions. We then turn to the TBA equations and the properties of their solutions. Finally we 
compare the results of the two methods and find the relation between UV and IR parameters. 

2 Boundary conformal field theory 

We begin with a sketch of boundary conformal field theory in general, although we shall often 
specialise to those theories described in the bulk by a minimal Virasoro theories with 'A'-type 
diagonal modular invariant. At the end we give the details for the Yang-Lee model. For a 
general discussion of the boundary content of the Virasoro minimal models see |j. 

A boundary conformal field theory is specified by giving the field content in the bulk, a 
list of the possible boundary conditions and the field contents of these boundary conditions, 
the fields which interpolate them, and some constants. 

In the bulk there are primary fields <& j of weights hi , hi which for a diagonal modular 
invariant are equal; under conformal transformations they obey 

*i(z,z) = 

where Xi = 2/ij. 

For the A-series, the conformal boundary conditions (which we shall denote by a) may 
themselves also be indexed by the Virasoro minimal representations (l|,|5[. We adopt the 
notation <p^-{x) for a field of weight hj which interpolates two boundary conditions a, /3, 
which may be the same or different. Such a field exists if the fusion coefficient N a p l is one, 
and does not exist if it is zeroQ We take such boundary primary fields to transform under 

'Remember that for A- modular invariant Virasoro minimal models all representations are self-conjugate 
and all the fusion coefficients are either 1 or 



dw 
dz 



$i(w,w) , 
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conformal transformations as^ 



dw 



dz 



d a P) ( \ 



Inside correlation functions we assume that the product of fields can be replaced by their 
operator product expansion - we do not try to interpret them as genuine expansions of the 
actions of operators, nor do we consider the associated problems of orderings of the arguments 
of the fields. 

There are three different operator products - bulk-bulk, bulk-boundary and boundary- 
boundary. The first are the standard expressions^ , 

$i(x,z) $j(w,w) ~ ^ Cij k \z - u,! 3 *-**- 3 * $ k (w,w) + . . . . (2.1) 
k 

Next we have the operator products of the boundary primary fields: 

k 

Finally, a primary bulk field near a boundary of type a may be expanded in boundary fields 

as 

$i(z,z) ~ Yl l 2 ( z " w )\ hj ~ Xi ^ Qa) H + • • • . (2-3) 

3 

where z is in the bulk and w is the point on the boundary closest to z. 

To determine the correlation functions consistently, one needs to normalise the different 
boundary conditions differently. One way to do this is to give the zero-point functions on a 
disc of radius 1 with boundary condition a, which we denote by Z a , 

Z a = (l)f iC . (2.4) 



cf>W\z) ( j>f' r \w)~J2 cia ij r)k \z-w\ hk ~ h ^ 4° 7) H + ... . (2.2) 



The rules for calculating the various constants in equations ( 27L|) , (p.2j), ( |2.3j ) and (2.4) are 



] , but some work is required 



essentially given in a series of papers by Cardy and Lewellen 
to obtain explicit expressions. For the structure constants of the A-type model, see ||. 

Using these structure constants, zero-point functions and the chiral blocks, we can assign a 
value to any n-point function, i.e. to any surface with some assignment of conformal boundary 
conditions to the boundaries, and some insertions of bulk and boundary fields; if the boundary 
conditions change at some point along a boundary, then there is of necessity a boundary 
changing field inserted at this point. 

We can now express various simple correlation functions on a disc of radius 1 in terms of 



t There is some choice in the transformations of boundary fields, see || 
'''n.b. in the 'A'-type invariant we only have spinless bulk primary fields 
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these structure constants: 



4 aa \e ie ) 



disc 
a 

disc 



disc 



^(re^)^p(l) 



(act) 
3 

t p) {6i) <j>f a \e 2 ) 



a 

disc 



disc 



o , (i )f c = z a 

{ { ^B}Z a ) (l-|*| 2 )~ x< , 

((*) B { c {a ™ h Z a ^ |l-2rcos0+r 2 |^ \\- r 2 \ h ^ Xi , 



(a/3a)l 



( A<*Pi)k ^( Q 7")i 
\ 1 



2sin%^ 



-2hi 



2sin^f 



° kk 
h k — hi — h 



z, 



2sin^f 



hi—hj—hk 



2sin^ 



(2.5) 

where 6ij = 6i — 9j, etc. 

We are interested in calculating correlation functions on the cylinder but we will be able 
to relate all the quantities of interest to simpler correlation functions on a disc. 

Let us consider the partition function of a right cylinder of length R, with the two ends 
being circles of circumference L with conformal boundary conditions a and f3. This can be 
realised as a rectangle in the upper half plane with vertices 0, L, iR and L + iR, with the two 
vertical sides identified. The partition function can be written as 



Z^(R,L) = Tr H(a0) (exp(-LflJ$(*))) . 



(2.6) 



By mapping this rectangle to a half-annulus in the upper half plane by z i— ► exp(7r z/R), we 
can use the operator formalism in the upper-half plane. The Hilbert space Ti.( a f3) decomposes 
into a direct sum of irreducible representations of one copy of the Virasoro algebra 



n 



and the Hamiltonian is given as 



(a/3) 



£>strip 



^ (L °"2l ) ' 



so that 



Na^xM , 



(2.7) 



(2.8) 



(2.9) 



where Xi{o) 1S * ne character of the representation i, q = exp(— tvL/R), and N a p l are the 
Verlinde fusion numbers. 

We can also map the cylinder into a complete annulus of inner and outer radii exp(— 2itR/L) 
and 1 respectively by z^w=exp(2i7rz/L), and use the operator formalism on the plane. The 
Hilbert space for the bulk theory on the plane carries a representation of two copies of the 
Virasoro algebra, in terms of which the cylinder Hamiltonian is given as 



- T (L 



(2.10) 
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This space carries a representation of two copies of the Virasoro algebra. The highest weight 
states | i ) are labelled by the highest weight representations of the Virasoro algebra and given 
by the action of the bulk primary fields on the S'L(2)-invariant vacuum | ) by 

\i)= lim <f>i(w,w)\0) , (2.11) 

w— >0 

and hence 

= %ai 1 (0|0) . 

We shall not necessarily normalise the vacuum to be norm 1, and indeed for the Yang-Lee 
model it is convenient to take ( | ) = —1. 

The partition function on the cylinder can be calculated in terms of boundary states | a ), 
| j3 ) in the Hilbert space for the bulk theory on the plane: 

Z% = {a\e- Rn ^\P). (2.12) 

For a boundary state to represent a circular boundary with conformally invariant boundary 
condition, it must satisfy 

[L m - L_ m ) | a ) = , ( a \ [L m - L_ m ) = . 

A basis of states satisfying these conditions are the Ishibashi states. There is one such state 
for each Virasoro representation in the minimal model in question, which we denote by 

1 + |0)) .| 1 + ^ + ...} |0) . 

«i| = (i|(i + ¥ 1 + ...} ■ «0| = <0|(l + ^# + ... 



2hi J " 1 N 1 t c/2 

These are normalised so that 

((i\e- RScilc M\j)) = 5 ijX i(q)(i\i) , 
where q = exp(— 4-nR/L). We expand the boundary states in terms of the Ishibashi states as 

I « ) = 9i I 3 )) ■ 

j 

(It is always possible to normalise the fields and states such that the coefficients g l a are real, 
and we shall assume that has been done. See Q for details). Using these states, the partition 
function is given from (|23| ) by 



Z (a/3)( r l' r 2) = 2^«^Xi(5)(»l*> ■ ( 2 - 13 ) 

i 

Noting that 

= (5). 
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where is the matrix implementing the modular transformation, we arrive at Cardy's result, 

gWf t W)=Y,n a js j i . (2.i4) 
j 

These equations can then be solved to find the g l a . 

One can also show that the same boundary states | a ) describe a circular boundary of 
radius one in the plane^Jand so the constants Z a can now be found by expressing the one-point 
functions on a disc of radius 1 in two different ways as 

(^Bf Z a = (M0))T c = (a\i) = Y J 9i((j\i) = 9i(i\i) , (2.15) 

j 

to find 

Z a ^Bf=gi(i\i) . (2.16) 

This completes the identification of the structure constants and zero-point functions. To end 
this section we recall how to find the conformal values of the ^-functions in terms of this 
conformal data. To fix the g- functions g a , we note that the leading term in the cylinder 
partition function is given by 

Z^faQ-fifle-^Wn) . (2.17) 

where | Q ) is the state of lowest conformal weight in the theory defined on the plane. This 
means that the ^-functions g a satisfy 

9a9p = 9a9?(^\^) ■ (2-18) 

However, recalling that the conformal vacuum representation '0' satisfies No a l = 5 l a , we see 
from ( gig ) that 

^ = o^)' 9 * 9 <* = Tnfnj- (2 - 19) 

Finally, we recall that the modular 5-matrix for the Virasoro minimal models satisfies Sn a >0 
for all representations a, and so if we require the (^-functions to be positive and the g l a to be 
real we should take ( £1 | O ) > and can read off their values from eqns. ( |2.18| ) and ( p. 19 ) as 

ID 1 1/2 Sfta n 9a /„ „m 

91 = lSno1 ' 9a = \^> 9a = 7wW ( } 

We now give the explicit results for the case of the Yang-Lee model. 



2.1 The Yang-Lee model 

The Yang-Lee model is the simplest non- unitary conformal field theory, il^s , and has central 
charge —22/5 and effective central charge 2/5. There are only two representations of the 
Virasoro algebra of interest, of weight and —1/5, and consequently only two bulk primary 

§ n.b. the normalisations of annulus partition functions differ from those of cylinder partition functions Jlo| , 
so some care must be taken to define boundary states for discs of general radius. 
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fields, the identity 1 of weight 0, and ip of weight = —2/5; equally there are only two 
conformally- invariant boundary conditions, which we denote by 1 and <£. The fusion rules 
are 

1x1 = 1 , 

1 x p = p , (2.21) 
ip x <p = 1 + <p . 

Prom these rules we read off that there are only three non-trivial boundary fields, all of weight 
he/, = 1/5. Two of these interpolate the two different conformal boundary conditions, 

which we shall both denote by tp, and one lives on the 3> boundary, 

<t> = ■ (2-23) 

As indicated^], we shall often simply denote the non-trivial boundary primary field on the 
boundary by <f>. 

If we label the two representations of the Virasoro algebra by (j) and 1 then the modular 
5-matrix is 

Su S 1(f) \ _ 2 f~ sin ir sin f \ ( -0.8506.. 0.5257.. 



S V% S WJ ^ V sinf sinf J V 0.5257.. 0.8506.. J' (2 ' 24) 
2.1.1 Yang Lee structure constants 

We give here the structure constants appearing in all the operator products of interest, that 
is the bulk OPE 

tp(z, z) ip(w, w) = C W A \z - w\ 4 / 5 + \z - w\ 2 ^ 5 tp(w,w) + . . . , 

the boundary OPEs, 

<Kz)<f>(w) = C^ 1 \z - w\ 2 / 5 + Cp/ \z - wl 1 /* <f>(w) + ... , 

4>{z)^{w) = C^z-w^ipH + ... , 

and the two bulk-boundary OPEs 

V (z)U = ^Bl\2{z- W )\^ + ... , 

¥>(*)!* = ^Bl\2 {z-w)\ 2 "> + ^B^z-w^^w) + ... . 

" In this model, there is only a single non-trivial bulk field, a single boundary condition admitting a relevant 
perturbation, and a single non-trivial boundary field, all labelled by the same Virasoro representation, and 
hence we have denoted them by tp, <I> and (j> respectively, in the hope that this will reduce confusion. 
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We want all these structure constants to be real, and a suitable choice is 





2 


1/2 


l+v^ 






1/2 


2 




2 


1/2 



Ca 



a 



a 



■ a 



a 



Mb* 



(#)£l 



r(i/5)r(6/5) 
r(3/5)r(4/5) 



1/2 



l+v^ 

1+v^ 
2 

5+^5 
2 



1/2 
3/2 
1/2 



(2.25) 



a 



The state of lowest conformal dimension in the bulk theory is 

\n) = \<p) , 

and we choose to normalise the bulk highest-weight states as 

<0|0> = -1 (<p\<p)=l. 



(2.26) 



(2.27) 



With this, and demanding gt, > 0, the coefficients appearing in the boundary states are 



9l 
9% 



V5+1 
2y/l 

VE-2 



1/4 
1/4 



9i 

9% 



■v/5-l 
2^5 

2+V5 
VE 



1/4 
1/4 



The zero-point functions are then given as 

Z 1 = (l)f c = (IL|0) = - 5 S, Z* = (l) 



disc 



($|0) 
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Finally, the (^-functions are given by <? a = \J ( <p \ <p ) = <7*, i-e. 

log 5i = 3 log 
log g<s> = \ log 

2.1.2 Yang Lee correlation functions 



\/5-l 
2Vf> 


= -I log 


l+v^ 

2 


2+VE 


= f lo g 


l+v^ 


VE 


2 



| log 5 
slog 5 



0.321482.. , 
0.159729.. . 



(2.28) 



(2.29) 



(2.30) 



In terms of the structure constants fl2.25| ), the simple correlation functions (2.5) on a disc of 
unit radius are 



/ -i \disc 
\ la 



■ disc 



, disc 



< 0(1) > 

<p(re«)#l) > C , 



disc 



Z a , ( 0(e 
C ffl )fl£ Z Q (1 - r 2 ) 2 / 5 
C^ 1 Z $ (2sin(^/2)) 2 / 5 







< ?(re") >™ 

dif 

C^ 1 Z$ (1 - 2r cos + r 2 ) 1 /5(l - r 2 )V 5 



(2.31) 



We shall also need two more correlation functions, which are given in terms of chiral four- 
point blocks. A basis for the chiral four-point functions of the weight —1/5 field are f 1 and 
f^, given in terms of hyper geometric functions by 



f\x) = [x{\ - x)] 2 / 5 F(l |; | ; x) , f*(x) = x 1 " (1 - x) 2 / 5 |; f ; x) . 



(2.32) 
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In terms of these chiral blocks we can give the bulk two-point functions on a disc of radius 1: 

(^(OMre^))" 1 " = Z a {c w VV) + (a) ^V/V)} , (2-33) 

which is rotationally invariant. We can also give the l-bulk-2-boundary point functions in 
terms of these two chiral blocks, but it is more convenient to perform some transformations 
on the arguments to make the expressions manifestly real, in which case we find 

/ \ disc c 

(<p(0) 0(1) # e »)) # = { (^(2sinf) 2 / 5 (cosf)- 4 / 5 F(A,| j; ii; _ t an 2 f) 

+ (^^/(2 S in|)V5 (cos | ) -3/5 jF( ^^.^. _ tan 2 f) } . 

(2.34) 

3 ^-functions and perturbed conformal field theory 

We define the (/-functions of the perturbed theory as 

(a\n) 

a ~ {Q\ny/ 2 ' [6 } 

where | a ) is the (possibly perturbed) boundary state and | Q ) is the ground state of the 
perturbed theory. Typically these functions will not be equal to the ^-functions g a but will 
differ by the extensive free-energy terms discussed in the introduction: 

\o g g a = logg a - Lf* cit . 

However it is the Q a that are accessible in conformal perturbation theory, and not the g a . 
Currently it is not possible to determine /a Cft directly from conformal perturbation theory, 
but using the truncated conformal space approximation (TCSA), it is possible to approximate 
the functions Q a numerically for a reasonable range of L, and this will enable us in some cases 
to estimate /« cft . For the case of the Yang-Lee model, the boundary conditions of interest 
are 11 (which admits no relevant perturbation) and &{h), the perturbation of the conformal 
$ boundary by the integral along the boundary 

h J <p(x) dx . (3.2) 

The bulk perturbation corresponds to a term in the action 

A J Lp(w,w) d 2 w . (3-3) 

The functions Q a therefore have expansions 



log a* (A, L) = £~ eUA£ 12/5 ) n , 



(3.4) 
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In section 3.2 we discuss the TCSA calculation of Q a , but before that, in the next section we 



calculate the following coefficients in their expansions: 

log£j(A,L) = log 5! + di (AL 12 / 5 ) + ... , 
lo g g m (X,L) = log 5 $ + c w (hL^>) + c 20 (/iL 6 / 5 ) 2 + c i(AL 12 / 5 ) + ... 

To the order to which we will work, the bulk and boundary perturbations can be treated 
independently, so we shall first calculate coi and d\, and then cio and C2o- 

3.1 The bulk perturbation 



For a purely bulk perturbation ( |3.3[ ), after mapping to the plane, the Hamiltonian of the 
perturbed theory is 

H circ (X, L) = H + A ill 

= ^(Lo+Io-^) + A jy(e ie )&9. (3.6) 

The full ground state is given by 

|0) = \^) + \Y^^a\^a)+0{\ 2 ) , (3.7) 

a 

where the sum is over an orthonormal eigenbasis of Hq (where the ' indicates that the ground 
state \ ip) is to be excluded) and the coefficients VL a are given by first order perturbation 
theory as 

n a = . (AlEihl . ( 3.8) 

(ip\H \(p) - (lpa\H \^ a ) 

Combining ([^) and ( |3.6[) , we find 

/ / T \ 12/5 1 f 

£ n a \ ^ a > = - - (l-P) - (l-P) j> ^e w ) \<p)M, (3.9) 

„ \ 27T J L + L + 2/5 J 



where P = \(p)(tp\ is the projector onto the state From ( |3.7D , the norm of | ) is 

unchanged to order A, and hence the correction to the (/-function is given by 

l ^7^l\ = A£ 12/ V + 0(A 2 ) , (3.10) 

where 

a a = -l(2vrr 12 / 5 / dg( a | (l-P) 1 (l-P) y (e") | y ) . (3.11) 
7 L + L Q + 2/5 

Using the rotational invariance of the boundary state and of ( (p | and the relations 
f 1 1 

x a dx = , x Lo+Lo cp(l)x- L °- Lo =x- 2/5 ip(x) , (Q|P = 5 $(yj|, 

Jo « + 1 
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we arrive at 



a a = -(2n)- 7 / 5 £ ^(±.{ a \<p(x)\<p) _ (^(x)^)) . (3.12) 
The correlation functions in ( 3.12| ) are 

(a\<p{x)\<p) = (^(xM0))f c = Z a {c^f\x 2 ) + ^i?} V/V)} 
(<p\<p(x)\<p) = -C^x 2 l\ 

so that using (HID, the expression for a a simplifies to 

a a = -{2^l\c^h + ^h), (3.13) 
where the integrals I\ and I2 are given (by Mathematica) as 

h = 0/H,V(dx^ ( l_ I 2 ) 2/5 F( 3 i 4.6. J 2 ) 



'0 x 

- r(i)r(| 

10sin(27r/5) T(|) 2 
and the somewhat obscure expression 



2.48171... , 



h = f^(/V)- 2/5 ) = £^((l-x 2 )^ ( |,|;f;, 2 )-l 



= -15,R,fl Z8-9 Ig-D + S ( 5 +v / 5) 3/2 -Vllogl+^I + 5n og5 _2) 
56 3-f2V- L ) 5)5)51 5 ' ' 8 v^tT 4 iu & 2 ^ 8 V lu & ° z / 

= -0.083937990... . (3.14) 
Substituting these values into a a we find 

di = oi = -(2tt)- 7 / 5 (C w ^/ 2 + ^/j) = -0.25311758.., (3.15) 

coi = a$ = -(2vr)- 7 / 5 (C w v / 2 + (^/l) = 0.0797648257... (3.16) 
While the expression for 1% is rather unwieldy, it cancels exactly from the ratio 

" log|-(A i ^)(2.)-'/ W (^_^) + ..., (3.17) 



so that CPT predicts 

c 01 - d x = (2tt)~ 2/5 5~ 3/4 - V5 ^ V5/ . - = 0.3328824.. . (3.18) 

' ' ' " " ' Y 

In section H we shall compare this with the TBA results from section 0. 



^-2/5 K -3/4 r (g) r (f) 

(i + V§)r(|) s 



n 



3.2 The boundary perturbation 



The boundary condition 9(h) corresponds to the perturbation of the 9 boundary by the 
addition to the action of the term 



SS = h 



/,a 



Z 0(Z 



(3.19) 



To calculate c\q and C20 we can assume that the bulk theory is unperturbed which means 
that we can replace the ground state | Q.) by \<p). The perturbed boundary state is 



{9(h)\ = {9\[1 - 5S + l(5S) 2 + 

and so the coefficients in the expansion 

( 9(h) \ip)=g$(l + c 10 hL 6 / 5 + 020 {hL^f + ... 

can be written in terms of the simple disc amplitudes 

( 27r )-6/s r ** 



(3.20) 



cio 



C20 



i(2tt) 
2 



-12/5 





2tt 2tt 



.„ \ disc 

(O v (o) ) $ d# , 



9i> 



{e id ) 0(e ie ') p(0) )~~ d6>d0' . 
/ $ 



disc 



(3.21) 







The correlation functions appearing in ( |3.2l|) are given in section 2.1.2 . The two-point function 
( (p(0)4>(e l8 ) }^ sc is rotationally invariant and so the integration in cio is trivial, giving 



cio 



-(27T)- 



1/5 ^BtpCfa 1 



(2tt)- 1 / 5 5 1 / 4 



r(|)r(f) 



2cos(7r/5)r(|) 2 



1/2 



0.99777.. 



(3.22) 



It is a little more complicated to find C20, which is given as 



C20 = (27r) 7 ^ 5 |c ( / ) 1L /3 + 
where the integrals J3 and ^4 are given by 



(*) B 
1 n <p 



'4 / , 



(3.23) 



h = r(2sinf) 2 / 5 (cos|)- 4 / 5 F(A,^ ; ii._ tan 2| )de 
J 

r (2sin f ) 1 /5 (cos |)-3/5 F( ^, £. _ tan 2 f ) d# 



(3.24) 



which do not seem to have simple closed expressions. However, by making use of the substi- 
tution v = cos 2 (#/2) and the explicit forms for the structure constants appearing in ( |3.23j ), 
one can show that the specific combination of J3 and I4 in (3.23) is given by 



1 .r 
C20 = {CIOY 



(3.25) 
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Taking the logarithm of ( |3,20 ), we find 



lo g g m (0,L) = log 5 $ + c w hL 6 / 5 + (c 20 - \(c w f) ihL^f + ... , (3.26) 
so that the coefficients c\q and C20 in (p-5|) are given by 



cio 



(2tt)- 1 / 5 5 1 / 4 



r(|)r( 6 ^ 



2cos(7r/5)r(f)2 



1/2 



0.99777.. , 



C20 = ( ^-|J (cio) 2 • (3.27) 



The coefficients cio and C20 have also been calculated using the free-field construction in [22], 
with the same results that we find. 

3.3 Evaluating log(g Q ) using the TCSA method 

The truncated conformal space approximation (TCSA) introduced in [|ll|] gives numerical 
approximations to the lowest eigenvalues of the Hamiltonian of a perturbed conformal field 
theory by restricting the Hamiltonian to a finite dimensional subspace, working out the matrix 
elements exactly and diagonalising the resulting matrix numerically. 

Originally developed for systems on a circle, we outlined in [|| how this method can also 
be applied to a system on a strip with conformal boundary conditions a, (3 on the two ends 
and perturbed simultaneously by a bulk field <3? of weight x$ and by boundary fields cj>i and 
(ft r of weights x\ and x r . The resulting Hamiltonian is 

H(R,X,hi,h r ) 

7T 

8=0 

The matrix elements of H can be easily calculated numerically and the Hamiltonian diago- 
nalised on spaces of the order of 100 states. 

The partition function on a cylinder of circumference L can be estimated numerically as 
the trace over the finite-dimensional truncated space. 

Z a{hi)i3{h r ){ L ,RiX) =Tt(exp(-LH(R,X,hi,hr))) . 

More details will be given in the next section where the method is applied to perturbations of 
the 1 and <£(0) boundary conditions, both by bulk perturbations, and (for the $(0) boundary 
condition) at the boundary. In the presence of the bulk perturbation, the mass scale is set by 
the mass M of the single massive particle of the scaling Yang-Lee model, which is related to 



A by [12,1 



M = kA 5//12 , k = 2 19 / 12 y^F = 2.642944 . . . , 

V 5 5 / 16 r(2/3)r(5/6) 

We present the results for the (1,1) and (l,<3?(/i)) boundary conditions in the next two 
sections. 
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3.3.1 The strip with (1,1) boundary conditions 

The (1,1) system has no relevant boundary perturbations and so the strip Hamiltonian is 
simply 

IT 

h(r,x) = ^l -^ + \\%\ 2 - x *J$(ex V (ie))de) . 

9=0 

The spectrum comprises scattering states of a single massive particle of mass M and no 
boundary bound states. In figure |l] we plot the 15 lowest eigenvalues of H/M against r = MR 
(as estimated from TCSA with 29 states). 

We estimate the partition function by taking the sum over all the TCSA eigenvalues 
at the given truncation level 

Z(r,0=^exp(-L^ TCSA ), 

i 

where I = ML. For any particular Z, we expect that the behaviour of Z(r,l) for small r is 
dominated by the conformal point, while at large r truncation errors will dominate. If we are 
lucky, there will be an intermediate range of r (the scaling region) in which Z(r,l) exhibits 
the scaling form 

log Z(r, I) a{l) r + b(l) + ... . (3.28) 

From equations ( |1.1[ ), ( |1.3|) and (|3.1|) , we expect that 

E cyl 

a(Z) = l5"' 6 (0=21ogSi(Z). 

To test for the behaviour ( |3.28| ) we plot log(dlog Z(r, l)/dr) to see if, for given Z, there is a 
region in r for which it is approximately constant. In figure [2] we plot log(dlog Z(r, l)/dr) 
against log(r//) for various fixed values of I calculated using the TCSA with 19 states and 
with 106 states. For small values of Z we see the typical behaviour expected - as r increases, 
log(d log Z(r, l)/dr) flattens out to a scaling region in which it is approximately constant 
before truncation effects take over. We can be fairly sure that this is indeed the scaling region 
by comparing the behaviour for different truncation levels, and, as expected, increasing the 
truncation level extends the scaling region further into the IR. On close examination we find 
for small Z < 1 that the scaling region is centred on log(r/Z) ~ 1.21 for 19 states and on 
log(r//) ~ 1.46 for 106 states. Hence we shall take the TCSA estimates of a{l) and b{l) to 
be those at log(r//) = 1.21 and log(r//) = 1.46 for 19 and 106 states respectively. For larger 
values of Z > 1 the effects of truncation start to interfere for logr ~ 2.75 and it is not clear 
whether we have actually reached the scaling region or not. 
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Figure [l]: The first 15 eigenvalues of H/M for 
the (1, 1) system plotted against r from TCSA 
with 29 states. 



Figure g: Plots of log(dlogZ(r, l)/dr) vs. 
log(r/l) for l—2 n , with (from left to right) 
n = 4, 3, 2, 1, 0, -1, -2, -3, -4. Results are for 
106 (solid) and for 19 (dashed) states. 



In figure || we plot a(l) as estimated in this way from the TCSA with 106 states and include 
a plot of Eq 1tc (1)/M for comparison. It is clear from the graph that it is possible to estimate 
a(l) effectively this way. 

In figure |I| we plot b{l) as estimated in the same way. The onset of the linear behaviour 
(|1.4j) expected for large / is clear in this plot. In figure || we plot db(l)/dl and give for 
comparison the exact TBA prediction from eqn. (|420| ). It is clear that while the TCSA 
cannot predict the correct value, it agrees with the TBA to our accuracy, and so in figure || 
we plot the combination 



£/r A +iog^ csA 



in 



f TBA 



showing its interpolation between the UV and IR values. Finally in table || we give the TCSA 
estimate of the coefficients £/( n )(0)/n! of the series expansion ( |3.4|) from truncations to finite 
dimensional spaces and the extrapolation of these coefficients to infinite level: 





63 states 


75 states 


TCSA 
90 states 


106 states 


oo states 


exact 


d 


-0.321615 


-0.321558 


-0.321518 


-0.321580 


-0.3215 


-0.32148269.. 


di 


-0.25230 


-0.25252 


-0.25269 


-0.25260 


-0.253 


-0.25311758.. 


d 2 


0.076715 


0.077208 


0.077523 


0.077570 


0.0775 




d 3 


-0.03510 


-0.03569 


-0.03599 


-0.03606 


-0.0360 




di 


0.01787 


0.01841 


0.01857 


0.01893 


0.019 





Table | 
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Figure |3[ a(l) as estimated from TCSA with Figure [|: 6(7) as estimated from TCSA with 

19 states (dashed line) and 106 states (solid 19 states (dashed line) and 106 states (solid 

line) together with Eq 1tc (1)/M (points) line) plotted against I. 
plotted against I. 



-o . 1 



-0.3 




Figure^ db/dl as estimated from TCSA 
with 19 states (dashed line) and 106 states 
(solid line) plotted against I. Also shown is 
the TBA prediction -2fJ BA . 



Figure^: log g%(l) vs. log? from /™ A and 
the (/-function from the TCSA with 19 
(dashed line) and 106 (solid line) states. Also 
shown are the exact UV and IR values. 



3.3.2 The strip with (l,<J>(/i)) boundary conditions: massive case 

We now consider the strip with boundary conditions (1, &(h)), corresponding to a Hamilto- 



man 



7T 

H(R,x,h) = ^(L -^ + \\§\ 12/5 J<p(ex P (ie))de + Hf| 6/5 <Ki; 



(3.29) 



The general partition function is now a function of three dimensionless combinations of R, 
L, A and h. In the massive case it is convenient to take these as the scaled lengths and the 
scaled boundary field, defined by 

r = MR, 1 = ML , h = h M~ 6/5 , (3.30) 

with the expectation that 

\ogZ {%Mh)) (R,L,X) ~ -ra(l,h) + b(l,h) , (3.31) 



b(l,h)=log(g l (L,X) 



^(% 3 +/i) + log(<7l(0s»$)(0)- (3-32) 
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To simplify matters, in this section we will only discuss the model with a purely bulk pertur- 
bation, so that h = 0. This corresponds to the value b = —1/2 of the IR boundary parameter 
defined in eqn. ( f4.4|) . In figure [?] we show the spectrum of the model with these boundary 
conditions and truncation to 81 states (the first excited state here is actually a boundary 
bound state). Examining dlog(Z(r, Z))/dr, we find that, as before, for small values of I there 
is a scaling region in r which broadens with increasing I before truncation effects intervene. 
For 140 states, the small / scaling region is centred on log(r//) ~ 1.38 while for large I trun- 
cation effects occur for logr ~ 1.5. For small I we can clearly see that scaling has set in and 
can be confident that we can estimate a(l) and 6(7) using the values at the midpoint of the 
scaling region. For large I we can assume that the 'best guesses' of a(l) and b(l) are given by 
their values at either (a) the small I effective value of r, or (b) the largest value of r before 
truncation effects are obvious. In other words we consider the two estimates of r e g(l) defined 
by 

(a) logr e g(l) = 1.38 + log / , and (b) logr e s(l) = min( 1.38 + log / , 1.5 ) . 

In figures || and ^ we plot a(Z) and d6/dZ using these two methods. We see that this second 
method of taking the 'best guess' at the large I values gives quite pleasing results for these 
two quantities while having the problem that the errors are hard to estimate. In figure 10 
we plot log(gi(/) <7$(o)(0) using method (b) and compare the results using the TBA value of 
f a and a 'fitted value' as was used in As we see, the large I behaviour of the ^-function 
appears to be improved by using a 'fitted' value of f a , but again the errors induced by doing 
this are unknown. 




Figure |7[ The first 15 eigenvalues of H/m for 
the (1, $(0)) system plotted against r from 
TCSA with 81 states. 



Figure g: a(l) vs. I, as estimated from TCSA 
with 140 states using methods (a) (dashed 
line) and (b) (solid line), together with 



E™ c {l)/M (points). 
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-0.02 - 
-0.04- 
-0.06 ■ 
-0.08 - 
-0.1- 
-0.12 - 




Figure 0: log(ffi(/)ff$( )(0) vs. log? from 
the TCSA with 140 states using method (b) 
and /™ A (dashed line) and a fitted value of 
f a (solid line). Also shown are the exact UV 
and IR values. 



Figure ® db/dl as estimated from TCSA 
with with 140 states using methods (a) 
(dashed line) and (b) (solid line), together 
with the TBA prediction. 



In table y we give the coefficients in the power series expansion ( |3.4| ) as found from the 
TCSA data by estimating the n-th derivatives of log(^i^$(o)) wr t A at A = 0. We also give 
the result of extrapolating these coefficients to infinite level. 

Combined with the results for d m in table |l] we can then estimate the coefficients (com — dm) 
which appear in the expansion of the ratio \og{Q^^/Qj). These are compared with the results 
from the TBA in section |5j, table [?]. 





81 states 


98 states 


TCSA 
117 states 


140 states 


oo states 


exact 


coo + do 


-0.16089 


-0.16108 


-0.16123 


-0.1613 


-0.162 


-0.161753565.. 


c i + di 


-0.1754 


-0.1750 


-0.1748 


-0.1746 


-0.174 


-0.173352755.. 


C02 + d 2 


0.06745 


0.06703 


0.06712 


0.06676 


0.0628 




C03 + d 3 


-0.0370 


-0.0370 


-0.0366 


-0.0363 


-0.033 




CQ4 + d4 


0.02248 


0.02267 


0.02230 


0.02225 


0.02 





Table 0: The coefficients (co n +d n ) from TCSA and CPT 



3.3.3 boundary conditions: massless case 

We now consider the (1, &(h)) model with the bulk massless. For this case, the three combi- 



nations (3.30) do not make much sense, and it is better to consider the system as a function 
of the two variables 

r = r/l = R/L, h = hl 6/5 = h L 6 / 5 . (3.33) 

As was shown in figure 1(b) of our previous paper Q, the spectrum of H here is real for h 
positive, and for h 3> 1 approaches that of the strip with (1, 1) boundary conditions. For 
sufficiently negative h, the spectrum becomes complex. 

The Hamiltonian on the circle is unperturbed if the bulk is massless, and so the ground 
state energy on the circle is simply given by ( 2.10| ): 



RE C Q 1TC 



TT R 7T _ 

(U) ~15 T ~ ~15 r ' 
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Similarly the massless limit of the linear term in Q, is simply given by (HH), so that 

( 3.31 ) and ( |3.32| ) become 



7T 



^gZ (1Mh)) (R,L,Q) ~ — r + b(h) , 



b(h) 



15 

2 -i/6 



h 



5/6 



+ log(#$(/ 1 )0M)) . 



(3.34) 
(3.35) 



In figures |Tl|a-p!l]c we show the TCSA estimates of Z(r, h) for various values of h for 
truncations to 29, 67 and 140 states. There is good convergence for r < 8. 





(a) h = -1 . (b) h = 1 . (c) h = 4 

Figure [Llf Graphs of log Z vs. r for fixed values of h from the TCSA to 29 
states (dashed line), 67 states (dotted line) and 140 states (solid line) 



In figures 12a-12~ we plot log(dlog(Z(r, h))/dr) vs. r to show the scaling region in which 
logZ grows approximately linearly, and include the exact value log(7r/15) from eqn. (3.34) 
for comparison. We see from these graphs that for small values of h there is good convergence 
for r < 3 — 6. 




-1.2 
-1.4 
-1 . 6" 
-1 . 8 



(a) h = -1 . (b) h = 1 . (c) /i = 2 . 

Figure [l2| Graphs of log(dlogZ/d?) vs. r for fixed values of h from the 
TCSA to 29 states (dashed line), 67 states (dotted line) and 140 states 
(solid line). Also shown is the exact value log(7r/15). 

We can estimate b(h) in two ways - using the exact value of REq 1vc or using the TCSA value 
calculated using eqn. (|3,32j). With b(h) calculated in either manner we can finally use the 



TBA value of (fL) to calculate the ^-function from eqn. (|3~35D . Infi gure |l3| we plot the TCSA 
estimate of log(g$(/i)) using the exact value of E™ c and the estimated value, together with 
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the massless TBA result (suitably normalised). We see that there is excellent agreement. 




-0.6 



Figure [l3| Graphs of log(g$(h)) vs. log(/i) from the TCSA with 140 states 
using the exact value of Eq" c (dashed line) and the estimated value (solid 
line) together with the TBA results (points). 

As a further check, in table ^ we give the TCSA estimates of the coefficients c n o- From [14|, 
we know that the TBA and CPT (/-functions agree in the massless case; in table | we give the 
values of cqq, cio and C20 from CPT and C30 from TBA. As can be seen, the extrapolations of 
coo and cio are really quite good but the TCSA under-estimates the value of C20 by about 15% 
and the value of C30 by a similar amount. However, the TCSA estimates are still increasing 
quite fast with truncation level and would probably agree better on increasing the truncation 
level. 





67 states 


81 states 


TCSA 
98 states 117 states 


140 states 


00 states 


CPT/TBA 


coo 


-0.161720 


-0.161731 


-0.161739 


-0.161745 


-0.161749 


-0.16176 


-0.1617535656.. 


ClO 


0.997714 


0.997731 


0.997742 


0.997751 


0.997758 


0.99777 


0.9977728216.. 


C20 


-0.044147 


-0.044428 


-0.044674 


-0.044890 


-0.045083 


-0.046 


-0.0525515368.. 


C30 


0.007897 


0.007965 


0.008024 


0.008076 


0.008122 


0.0084 


0.00966359* 



Table |: The coefficients c n0 from TCSA, CPT and TBA (denoted by *) 

4 The TBA approach to ^-functions 

The TBA is an alternative to conformal perturbation theory and is characterised by its use 
of IR scattering data - the bulk scattering matrix of the massive particles, and the reflection 
factors describing the scattering of the particles off the boundaries. 

Exact TBA equations describing the g-function flow for both massive and massless bulk 
were proposed in Q , following on from earlier work in [|l5|] . The derivation of these equations 



(which we recall briefly in section 4.2) involved a number of assumptions, but to date only 
some simple limiting behaviours have been checked. We will perform a more detailed analysis 
for the Yang-Lee case, which will lead to the conclusion that the equations of [||] do not tell 
the whole story in every situation. 



4.1 The Yang— Lee scattering data 

As a massive scattering theory, the Yang-Lee model can be thought of as a reduction of 
the Sine-Gordon model at coupling constant 8ir/^/ = 3/2 (in the normalisation of ju!]). The 
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bulk spectrum reduces to a single particle (the first Sine-Gordon breather) with two-particle 
S-matrix JTfl] 

In |||, the boundary reflection factors corresponding to the boundary conditions &(h) and 1 
discussed in the previous sections were identified as 

= R b (9) , Rj(6) = Ro(9) , (4.2) 

where 

MO) = (|) (|) (I)" 1 (S(9 + «r*F)S(* - ztt^))- 1 (4.3) 

These reflection factors are special cases of the general Sine-Gordon breather reflection factors 
proposed in |l8| with the parameters rj, i? of |l8| , |l9|| taking the values 

2rj ,b , 2i$ ,b 

V=± 9+ 2 > — = ± (o+ 1 )- 
7T Z 7T Z 

In j| the relation between b and /i was conjectured on the basis of numerical fits to be 

h = - \h c \ sin(vr(6 + .5)/5) M 6/5 , h c = -0.685289983(9) . (4.4) 

An exact formula for h c was also proposed, though no derivation was given; that gap will 
be filled in section 4.3.5 below. Taking b real gives |/i|<|/i c |M 6//5 ; to obtain h< — |/i c |M 6//5 or 
h>\h c \M e / 5 we take b=—3—ib and b=2+ib (with b real) respectively. 

The massless limit is obtained by taking ML = I — ► while keeping h constant, where 

h = hL 6 ^ 5 = - \h c \ sin(f (b + ±)) / 6/5 . (4.5) 

Since \h/{h c M)\ — * oo in this limit, b takes the values —3 — ib and 2 + ib for /i > and h < 
respectively. 

For /i > 0, with b = — 3 — ib, we have sin(^(b + |)) = — cosh(7r6/5), and for 6^0, 

K ~ l\h c \ (MLe 7 ^ 6 ) . (4.6) 

Hence the massless limit is I — > 0, 6 — > +oo, (/e 71 "^ 6 ) finite. Substituting 6 = — 3 — into ( |4.3[ ) 
we find 

*-3-*W = (I) (1) (ir 1 + t)s(e - f))' 1 , (4.7) 

which is still a pure phase for real 9. Hence, although b is now complex, Rb should still 
describe a physical reflection process. 

For h < 0, b = 2 + i&, so sin(f (6 + |)) = cosh(7r6/5) and for 6 » 0, 

K ~ -i \h c \ (MLe* % f*\ . (4.8) 

Hence the massless limit with h < is I — > 0, 6 — > +oo, (/e 71 " 6 / 6 ) finite. Substituting b = 2 + ib 
into (|PD , 

W) = (I) (I) (I)" 1 (SO? " f + + f - if )V 1 , (4.9) 



21 



which is not in general a pure phase for real 9. Hence, i?& should not describe a physical 
reflection process and (in particular) we can no longer expect to find a real spectrum in the 
IR limit, which was indeed what was found in [|| in the TCSA analysis of the spectrum in 
the massless case with h < 0. 



4.2 The TBA for ^-functions 

The starting-point for the TBA analysis of |3],[l5| was the generalisation of the Bethe-ansatz 
quantisation equation for particles on a line segment of length R — * oo. Far from the end 
points, the scattering is characterised by a two-body S-matrix element; information about 
the reflection of the particles at the end points is encoded in the reflection factors R a (9) and 
Rp{9). The main difference with the standard TBA analysis is that the quantity of interest 
is a sub-leading term (order 1/R with respect to the bulk quantities) of the free energy. 

The initial assumption is that a particle state can be approximated in the thermodynamic 
limit (R — ► oo) by a distribution of particles in scattering states of given rapidities satisfying 
self-consistency equations which simultaneously define a density of possible particle states 
p(9) and a density of occupied particle states p r (9), related by 

/oo 
p r (9')<l){e-e')&o' , (4.io) 
-oo 

m = -^ lo S^) > MO) = ~i^logR a (0) ■ 

Assuming that the number of allowed states with effective particle density p r (9) is M[p r (8)], 
one can write the partition function as an integral 

Z = e ~ LE = M/W] M[p r (9)] e - LE ^ pr ^ , (4.11) 

states 

where the energy of the configuration is E = J^°(M cosh 6) p r (0) d9 . The standard TBA 
method is to calculate N[p r {6)} by taking the number of allowed configurations in the interval 
A# with a fixed number of occupied states Rp r (9)A6 to be 

(Rp(e)Aoy. (412) 



{Rp r (9)A6y.{Rp(9)A8 - Rp r (9)A6)\ 



and replacing the factorials by the two leading terms logT(z) ~ zlogz — z + . . . of Stirling's 
formula, so that the total number of configurations with a given effective density p r (9) is 

AV(0)]~exp / (p\ogp-p r logp r -(p-p r )log(p-p r ))d9 . (4.13) 
L Jo J 

We can then calculate Z in the limit R — > oo by the saddle point method, giving the leading 
behaviour 

log Z ~ -RE^(L) + log(g a (L) g p (L)) , 
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where (extending the range of 9 by symmetry where necessary and setting e = \og{p/ p r — 1) ) 

/oo 
mcosh#L(#)^ 

" , (4-14) 

\og{g a (L)g p {L)) = / U a (9) + <j> p (9) - 20(29) - 2*6(6)) L(9) f 

J— oo ^ ' 

where = log(l + e~ e W) and e(0) solves the equation 

/oo 
0(9 - 9') L{9') g . (4.15) 
-oo 



This gives exactly the same result for Eq tc (L) as the usual TBA on a circle [12]. The g- 
functions are then identified as 



logGfo(L)) = fZ(M0) - 4>(20) - nS(e))L(9)§ 



(4.16) 



There are a number of possible sources of error in this expression. Firstly there are possible 
corrections from the next-to-leading term (— logz/2) in Stirling's formula, and secondly from 
the corrections to the saddle point evaluation of the integral in ( [4.11 ). Finally, the initial 
assumption that the Hilbert space for the model on a line segment can be satisfactorily 
approximated by scattering states might be wrong. Unfortunately we do not as yet know 
how account for these effects in a consistent way. In [|j they were formally bypassed by 
assuming that they give an overall contribution (also appearing in the periodic case) which 
can be neglected when considering changes in the ^-functions. Further problems may arise due 
to the presence of boundary bound states; again, these were neglected in the above discussion. 
In section 4.3. 1| below the corrections that are necessary for the Yang-Lee model in this case 
will be found empirically, but a first-principles treatment is still lacking. With these caveats, 
equations (|4.15|) and (f4,16|) form the prediction of |||. 



4.3 Properties of the TBA g-functions 



From ( 4.15| ) and ( 4.16| ) it is possible to derive several analytic properties of the TBA g- 
f unctions. 



In section 4.3.1 we examine possible ambiguities in the g-functions, and give the prescrip- 



tion which we have adopted in this paper. In section [4.3.2| , under the assumption that the 
TBA g-functions are the sum of a term linear in L and a term which is regular in L 6 / 5 , the 
coefficient of the linear term is evaluated analytically. This enables us to define the TBA 
(/-functions, which are directly comparable with the TCSA results. 

The behaviour of L(9) in the UV and IR limits is well understood so it is possible to find 
the net change in the (/-functions along the flow. The g- functions in eqn. ( 4.16| ) suffer from 
ambiguities due to the possibility to change the integration contour, but we show in section 



4.3.3| that there are no contours for which the net change in the TBA g-functions agrees with 
the TCSA res ult. 

In section 4.3.5| we show that the expression for the ratio log^^m) (0/51 (0) is rather 
simpler than that for the two (/-functions separately, as was also the case in the perturbed 
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conformal field theory analysis. This enables us to make an expansion in the massless case 
and obtain the analytic value of h c defined in eqn. (4.4). 

We give some numerical power series fits to the TBA results in section 4.3.6 , and find a 
combination of functions which appears to be regular in I 12 / 5 . 

Finally in section 4.3.7 we discuss the singularities that occur in <jwm for sufficiently 

negative h and relate these to the singularity in the ground-state energy on the strip observed 
in a. 



4.3.1 The contour prescription for the (/-functions 

If we examine the formula ([4.16| ) for the ^-function in a little more detail, we notice a couple 
of immediate problems. Firstly, the two distinct boundary conditions <£(/i(6 = 0)) and 1 are 
described by the same reflection factor and hence have the same kernel function fa, yet have 
different (/-functions. Secondly, 4>b describing the $>(h(b)) boundary condition has poles which 
are 6-dependent and (in particular) cross the real axis for b = ±1. This means that taking the 
naive interpretation of ( [4,16] ), with the integration contour always along the real axis, will lead 
to (/-functions which are discontinuous at b = ±1 whereas the physical (/-functions should be 
continuous in b. The resolution of these two problems appears to be that the contour in (|4.16[) 
depends on the boundary condition, and in particular is different for the 1 and &(h(b = 0)) 
boundary conditions. Changing the integration contour will mean encircling some of the poles 
of 4>b, the poles so encircled being called 'active', in analogy with the terminology adopted 



in [20 in a related context. 



The exact choice is, at this stage, a matter of trial and error. In section 4.3.3| we show 



that there are no choices which will give the same (/-functions as the TCSA, but with the 
following prescription at least the ratios of (/-functions agree: 

1. The 1 boundary condition 

We take \og[g%) to be given by the straightforward application of ( [4.16 ), that is 



log(si) = (<M#) - <K20) - it6{0j)L{6) 



de 

4tt ' 



with the integration contour along the real axis. 
2. The <3?(/i(6)) boundary condition 



We take \og(g^h(b))) to be given by (|4.16|) for — 3 < b < —1, that is 



log(<7$ (M &))) = f°° ( MO) - 4>{20) ~ ) L(0) g , -3 < b < -1 , 

with the integration along the real axis. For all other values of b, it is given by the 
analytic continuation of this function. In particular, for h>\h c \ (i.e. b=—3—ib) in the 
massive case, and for h > in the massless, the integration contour may again be taken 
along the real axis. 
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The simplest way to derive the continuation of log(g^^i l (b))) is to isolate the terms in the 
integrand which have 6-dependent poles and treat these separately. We have 



MO) = M0) - Ho - - HO + wr*^) , 

so that for —3 < b <— 1, 

log(3 W)) ) = log( M )- I" U{e-m b ^)+<t>(e + m b ^j)L{e)f. (4.17) 

When the integration contour is along the real axis (but not necessarily otherwise) the two 
terms in the integral ( 4.17| ) give the same contribution, so at least for — 3< b<— 1 we have 



/oo 
0(0_ i7r k£3)£(0)d£ (418) 
-oo 

The term in the integral (p{9 — iir^^-) has poles at 

9 = + 2nm , 9 = iir^± + 2nm , 9 = iir^ + 2niri , 9 = in^f- + 2nm . (4.19) 

This means that as 9 passes from b < — lto6>— 1, the contour must deform away from the 
real axis to avoid a discontinuity from the pole at 9 = iir(b+l)/6, or equivalently the pole at 
*7r(&+l)/6 is 'active' for b > — 1, and similarly the pole at m{b— 1)/6 is 'active' for b > 1, as 



shown in figure 14. 



Note that the poles that are 'active' are exactly those associated with boundary bound 
states 0,^]- It is clear that the transition of a scattering state into a boundary bound state 
is always associated with a pole in (f> a crossing the real axis, and in this case the converse is 
also true. (Recall that the general interpretation of simple poles can be subtle, and not all 
correspond to boundary bound states - see |pl| .) Hence, this seems the natural prescription 
for the integration contour in the presence of such bound states and it would be interesting 
to arrive at this result from first principles and check it in further models. 




-3 < b < -1 -1 <b< 1 l<b<2 



Figure |l4|: The determination of the contour through the analytic 
continuation of logQjw/^j,)')) from b<— 1 to b>l. 
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4.3.2 Analytic structure of the TBA (/-functions. 

We expect the (^-functions as evaluated in the TBA to be regular functions of L 6//5 after the 
removal of an 'extensive' part, linear in L. Thus they should have the form 

log g a = Lf a + log Q a . 

Under this assumption, the coefficient of the linear term can be evaluated analytically as 



fa 



M-g a {l) 



1=0 



Al. Zamolodchikov showed in |12j that, in the limit / — > 0, > 0, 



L(0)~Lki nk (0-logf) 



with 



/_ 



oo 



de 



V3 



Substituting this into (|4.15|) we arrive at 

/ TBA = ^ +gin ^ M 



(4.20) 



independent of the contour prescription we choose. In the massless limit, using Q4.6| ) we find 
for h > that 



Jb L 



2 -l/6 



h/h c 



5/6 



(4.21) 



4.3.3 Net change in the TBA ^-functions along the flows 

Given (gig ), it is straightforward to calculate log(g a ) in the two limits L — > and L 
using the result 

log(i^) , UV, 
, IR. 

While the IR result is simple, 



L(9) 



(4.22) 



log ( 0a ) 



IR 



(4.23) 



the UV limit depends on the choice of integration contour in ( |4.16| ) . Although we have given 
above a particular prescription for the contours, no matter which contour is chosen, one must 
end up with 



log (5a ) 



uv 



flog 



n E 



Therefore the most general prediction of ( [4,14| ) would be an interpolating flow 



log (g a ) 



uv 



log (g a ) 



IR 



log 



n E 



(4.24) 



(4.25) 



There are two distinct physical situations: massless bulk and massive bulk. 
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1. Massless bulk 

There is a single identifiable massless flow with h > 0, from $ to 1, for which (|2~30[ ) 
gives 

logGfc) - log(si) = log(^) , (4.26) 



and which is in agreement with ( |4,25| ) for n = 4. 

2. Massive bulk 

For the massive flows from 1 and from with h>- 



\h c \, guv are given by conformal 
field theory, and our expectation that for the massive theory g\iR = 1 (already proposed 
in [Q] and well-supported by the TCSA results) leads to the results 



log (5i) 
log (5$ (ft) ) 



uv 



uv 



log(s) 
log (g) 



IK 



IR 



4 log 



flog 



1+V5 
2 

2 



I log 5 , 
I log 5 . 



(4.27) 



These are in contradiction with ( 4.25| ). 



Thus we see that the TBA equations ( |4. 15 ) do not agree with the TCSA even at the level of 
the overall change in the ^-functions without investigating finer detail. 

We note however that the ratio log(g<j>(/i)/<?i) also changes by log(^(\/5 + 1)) along a 
massive flow, so the TBA equations do have the possibility to describe the ratio correctly, 
and indeed we shall find supporting evidence for this later on. 



4.3.4 The relation of the TBA (^-functions to the Y function 

Having written the ratio of the g- functions in the form (4.18), we observe that the integral 
on the right hand side is exactly of the form appearing in the TBA equation determining the 
pseudo-energy e: 



log( 



/*(A(6))> 



log (51) 



ITT 



e(m^)+/sin(f 



(4.28) 



If we absorb the linear term into the ^-functions following ( f4.20D we obtain exactly the Q- 
functions, so that the ^-functions which follow from the TBA equations of [|| predict the 
particularly simple form for the ratio of the (/-functions 



log [G$(h{b)) I £1. 



l og y (ivr fe±3; 



(4.29) 



The same result was also obtained in by a different method, but with the restriction 
to the massless case. Given this, the extension Q4.29 ) to the massive case is perhaps not so 
surprising, but the derivation via the TBA seems unrelated to the methods of []l4j]. Just to 
finish this section, we comment that we have subjected (4.29) to several tests outlined in later 
sections, all of which support this identity. 
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4.3.5 The expansion of \og(g^ h Q ) ^(l) ) in the massless limit 

The perturbative expansion of e(6) in the massless limit has been found in [14 , 22]. If e{9) 
is normalised so that for large (lexp9) it has the asymptotic behaviour e{9) ~ 
small (I exp 9) it has an expansion 

e(*) = logi±^ + C.ile^ + C 2 (/e e ) 12/5 + ••• , 



then for 



where (from appendix A of |22|j ) 

2 -7/5^7/5 r( l) 



Cl 



Til) 



6/5 



° 2 ~ '75 ~5 



cosfr(|)r(f)2r(| 

From perturbed conformal field theory the coefficients in the expansion 



Q 2 - 



arc 



log- 



Si(0,L) 

and were found in section |3.2| to be given by ( |3,27| ) 



log^ + c 10 (hL^) + c 20 (hL^f + ... 
91 



Cl0 = (2 7 r)- 1 /5 5 V4 



r(!)r(f; 



2cOs(7T/5)r(|)2 

Given that in the massless limit we have 



1/2 



C20 



log- 



£i(0,L) 



O 



we deduce that 
ft fi = = 



-vr 3 / 5 2 4 / 5 5 1 / 4 - 



sin 



2- 



r(l; 



6/5 



(r(|)r(|))V2 

This is in perfect agreement with the numerical result obtained in 



0.68528998399118 . 



(4.30) 



(4.31) 



(4.32) 



(4.33) 



4.3.6 Series expansions of ingredients of the TBA expressions 

In an attempt to understand better the discrepancies between the TBA results and those 
from other approaches, we have investigated numerically the various ingredients of the TBA 
(/-functions. The aim was to determine their analytic structure by various power series fits in 
the variables x = AL 12 / 5 and y = /iL 6//5 . These results give highly accurate determinations 
of the first few coefficients, but the accuracy rapidly tails off. 

In each case we initially assumed an analytic expansion in x 1 ! 2 with the addition of 
possible non-analytic terms of the form a; 5 ™/ 12 , with the results given in table 0. On the 
basis of these results we then performed a fit to more restricted functional forms (with fewer 
non-analytic terms) with the results as in table ^. Finally, we know on theoretical grounds 
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that e(in(b + 3)/6) has an expansion in x and h and we give the numerical determination of 
these coefficients in table |(| 



-±L(0) = §e(0) - e(wr/3) 


-0.240605912529803- 2.2 lO-^x 5 / 12 + 0.419185703387 x 1 ' 2 
-3. 10" 12 a; 10 / 12 - 0.26930751 x - 3.6 10~ 8 a; 15 / 12 + 0.0120551 x 3 / 2 
-1. 10~ 7 a; 20 / 12 + 0.07 1 38 a; 2 - 2.5 lO^a; 25 / 12 - 0.0096 a; 5 / 2 + ... 




0.481211825059603-4. 10~ 14 a; 5 / 12 - 0.67825671567a; 1 / 2 

-9.1 10~ 13 a; 10 / 12 + 0.308598962 x- 7. 10~ 8 a; 15 / 12 + 0.05575394 x 3 / 2 

-1.3 10~ 7 a; 20 / 12 - 0.08 1 08 a; 2 - 2.6 10~ 5 a; 25 / 12 - 0.023 a; 5 / 2 + . . . 


-^{26) L{6) ^ 


0.240605912529 - 7.3 10~ 10 a; 5 / 12 - 0.67825671 x 1 ' 2 
+0.8102843 a; 10 / 12 - 0.33287 x - 0.0002 a; 15 / 12 
+0.062 x 3 ' 2 - 0.045 a; 20 / 12 + . . . 




-0.4812118250596 + 0.967384888 a; 5 / 12 - 4.3 10~ 9 x 1 ' 2 
-0.8102843 a; 10 / 12 + 0.20572a; + 0.00026a; 15 / 12 
+ . . . 


/_"(*)(9)-^9p(9)i 


-0.2406059125298+ 0.96738488748 a; 5 / 12 - 0.6782567157 x 1 / 2 
+6. 10" 9 a; 10 / 12 - 0.1271498 x + 5. 10~ 7 a; 15 / 12 + 0.055748 x 3 ' 2 
+9. 10~ 6 a; 20 / 12 - 0.02 3 50 a; 2 + 0.00004 a; 25 / 12 + . . . 



Tabic H the results of inital fits to the TBA data 



-§L(0) 


-0.240605912529802+ 0.4191857033888 a; 1 / 2 - 0.269307520 a; 
+0.0120552 x 3 ' 2 + 0.07138 a; 2 - 0.0096 a; 5 / 2 
+ ... 




0.481211825059603- 0.678256715680 a; 1 / 2 + 0.30859896278 x 
+0.055753953 a; 3 / 2 - 0.081082 a; 2 - 0.023 a; 5 / 2 
+ ... 




0.240605912529805 - 0.67825671 x 1 ' 2 + 0.8102848 a; 10 / 12 
-0.332875 x - 0.00020 a; 15 / 12 + 0.0615 x 3 ' 2 
-0.043 a; 20 / 12 + ... 




-0.481211825059593 + 0.96738488736 a; 5 / 12 - 0.8102849 a; 10 / 12 
+0.20572 x + 0.00021 a; 15 / 12 - 0.0060 x 3 ' 2 + 0.043 a; 20 / 12 
+ ... 


Oto(<?)-#20))L(0)g 


-0.240605912529802 + 0.9673848874698 a; 5 / 12 
-0.6782567156814 x 1 ' 2 - 0.127149765 x + 0.05575395 x 3 ' 2 
+0.023491 a; 2 - 0.020 a; 5 / 2 + . . . 



Table 0: the results of more restricted fits to the TBA data 



29 



e(m±P) = mLsin f - f~ x ^in b -^ - 0)L{0) d9/(2n) 



0.48121182508 + 0.332882394 x - 0.09242282 x 2 + 0.040463 x 3 - 0.02047 x 4 + 0.0104 x 5 - 0.00385 x e 
-0.9977728224 y - 0.086482609 xy + 0.03812851 x 2 y - 0.01991664 x 3 y + 0.010374 x 4 y - 0.00377 x 5 y 
-0.052551643 y 2 + 0.02553410 xy 2 - 0.0148490 x 2 y 2 + 0.008372 x 3 y 2 - 0.003206 x 4 y 2 
-0.00966359 y 3 - 0.0076724 xy 3 + 0.005150 x 2 y 3 - 0.00225 x 3 y 3 
-0.0021742 y 4 + 0.002199 xy 4 - 0.00119 x 2 y 4 
-0.000563 y 5 - 0.00046 xy 5 - 0.00016 y 6 



Table ^ numerical fit to the expansion of e 



We immediately notice that with the contour prescription of section 4.3.1, not only does 
<7i not have the correct net overall change, but neither is Q\ a pure power series in x (A). 
Even sacrificing universality, with the lack of such a basic property there is no hope to identify 
this quantity with a properly defined g-function. 

However, we observe that it is possible to find a function which has a regular expansion in x 
from the various ingredients of the TBA. If we drop the term —L(Q) /2 (which was introduced 
in H to account for the absence of a zero-momentum single particle state on the strip), and 
further, instead of running the contour on the real axis, we instead also encircle the singularity 
of <pQ at 8 = in/ 6 clockwise to pick up an extra contribution —L{m/Q) = —e(iir/2), we arrive 
at the function 

for which we find the expansion 



(0o(0)-0(20))L(0) — -e(f), 



-0.7218177375894 + 0.967384887468 x 5 / 12 + 2. 10~ 13 x 1 ' 2 - 0.435748727 x 
+ 6. 10~ 10 x 3 / 2 + 0.10457 x 2 + 0.004x 5 / 2 

which has (within our accuracy) better prospects of being a purely analytic function of x, 
plus a term in x 5//12 which we can identify as (\/3 — 1)ki 5 / 12 /2 = (\/3 — 1)1/2. 

To conclude, the g- functions predicted by the TBA equations of Q do not match the 
quantities derived perturbatively in this paper, nor do they have the expected analytic prop- 
erties. Although we have partial evidence that the latter problem could be overcome with an 
ad-hoc contour prescription, at the present stage of our understanding we have no definite 
justification for the above-mentioned discrepancies. 

4.3.7 Singularities for h < —\h c \ 



In section 4.3.1, we made the following identification in the case of massive bulk: 

log(<7*(h(6))) -logfai) = e(m b -^) +/sin(f ) . (4.34) 

The TCSA analysis in 0] revealed the existence of a critical value of /iM -6 / 5 , in the range 
—0-8 to —0-6, for which \og{g^n l n ) y\g%) diverges — this was signalled by a pole in a rational 
fit to the ^-function. Furthermore, the ground state energy on the strip has a square-root 
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singularity at some real value of r for /iM~ 6 / 5 < h c M~ 6 / 5 = -0.68529.. (see section 6 of Hi). 
This led us to the conjecture that the singularities in the two quantities occur at the same 
critical threshold /i c M -6 / 5 . Armed with relation ( [4.34 ) we can now prove this. First note 
that the quantity Y(9,r) = exp(e(#,r)) is entire in 9; further, for r = RM real, it has zeroes 
on the line Im.9 = 5ir/6 (or equivalently Im9 = — 5tt/6), arranged symmetrically about the 
imaginary axis. While their exact locations 9 n are not known, for large r a good approximation 
is sinh(Re# n ) = ±7r(2n + l)/r (n = 0, 1,2, . . .). Hence as r grows these zeroes move towards 
Re9 = 0. 

Suppose now that Re 6 = 2 so that the value 9(b) = in^^- relevant for ( [4 .34] ) lies on 
the line Im^ = 5n/6. As r increases, the zeroes 9 n of Y move along the line Im9 = 5tt/6 
and at some critical value r c (b) one will pass through 9(b), so that Y(9(b),r c (b)) = 0. This 
confirms the presence of logarithmic singularities in (|4.34[) , fixes the critical threshold at 6=2 
and proves (through eq. ( |4.4D ) the above mentioned conjecture. A similar analysis shows that 
logarithmic singularities appear in the massless h < case too. 



5 Comparison of TBA, CPT and TCSA results 



Firstly, it is a consequence of [14| that the TBA and CPT methods give the same answer for 
the (/-function along the massless flow. This is borne out by the comparison of our TCSA 
data with the TBA. However, we find disagreement between the TCSA and TBA predictions 
for the massive flows. 

On general grounds we showed in section |3| that the perturbative parts of the ^-functions, 
denoted by (/-functions, should have a series expansion for small L in L 6//5 . However, when we 



try to combine the power series expansions for the various ingredients given in section 4.3.6 
we see that the (/-functions predicted by the massive TBA equations of (3| do not have the 
correct behaviour. This is compounded by the fact that in the massive case the net change 
in the ^-functions does not agree with the CFT predictions. This leads us to our conclusion 
that the (7-functions and (/-functions given by the massive TBA equations of || do not agree 
with those obtained viewing the theory as a perturbed conformal field theory. 

One of the most natural sources of this discrepancy is from a neglected sub-leading term in 
the saddle point evaluation of log Z, but we have not been able to calculate such a term. What 
is clear is that the discrepancy between the TBA and CPT calculations of the ^-functions is 
independent of the boundary conditions. This is consistent with our observations in tables ^ 
and U that the first few coefficients in the power series expansions obtained from the TBA and 
TCSA for the ratio of (/-functions, log(Q^rf l \/Q%) do agree and that the extensive free-energy 
terms (the linear behaviour of the (/-functions for large L) are also the same in the TBA and 
TCSA results. 





TCSA (00 states) 


TBA 


exact 


coi 


-di 


0.332 


0.332882394 


0.332882407.. 


C02 


-d 2 


-0.092 


-0.09242282 




C03 


-d 3 


0.039 


0.040463 




C03 


-d 3 


-0.02 


-0.02047 





Table @: The coefficients (c 0n -d n ) from TCSA, TBA and CPT 
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As further evidence we show in figure 15 the TCSA and TBA estimates for the ratio of 
5&=-i.l and g b =-i.2- 




Figure ^ The difference log(p b= _i,i) - log(g 6= _i. 2 ) 
vs. log Z for the massive YL model. The solid line is the 
TBA result and the points are from the TCSA to 117 
states. 

6 Conclusions 

We have investigated the relation between the (/-functions proposed on the basis of a TBA 
analysis in || and those obtained from perturbed conformal field theory. 

We found that the massless case is correctly described by the TBA (/-functions (up to 
an overall constant) but that the massive case is not. However, in all cases, the ratios of Q- 
functions for different boundary conditions of the same model are correctly predicted by the 
TBA. A consequence of the TBA proposal of || is that this ratio of perturbative (/-functions 
should satisfy 



(6.1) 



This relation is already contained in the results of JhJ] in the massless case, but the derivation 
there seems rather different to ours. We have tested ( |6.1D in the massive case in two ways: 
we have compared the power series expansions from fits to the TBA data to the perturbed 
conformal field theory and TCSA calculations of these expressions; we have further compared 
the ratio of (/-functions directly over a range of I using the TBA and TCSA. Both tests 
support (|6T] ). 

We have also found strong support from the TCSA for the correctness of the TBA calcu- 
lation of the 'extensive' boundary free energy f a . 

To summarise, we have found strong evidence for the relation 



fJ, 



CFT (i) = 9 T a BA (D + Hi) 



where k(l) is an (unknown) function of I vanishing at / = oo, and which is independent of the 
boundary conditions, and where g™ A is the function which results from the equations of Q, 
analytically continued from the massless regime. In section |4.3.6| we further conjectured that 



„CFT 



(0 = 9, 



li!A (0 + |L(o)-i e (f) + MO, 
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where k(l) is an (unknown) regular function of Z 12 / 5 vanishing at I = oo, and which is in- 
dependent of the boundary conditions, and where again gJ BA is analytically continued from 
the massless regime. This is just the sort of relation that would arise if the TBA prediction 
of |H has neglected some corrections arising as next-to-leading order terms in the saddle-point 
evaluation of the path integral giving the partition function. As yet we have not been able to 
check this possibility. 

One of the unresolved problems we encountered was the correct contour prescription when 
the kernel function <j) a contained parameter-dependent poles. The result we found in section 
[4.3.1 is that it is consistent to take the contour to encircle those poles which correspond to 



boundary bound states. As mentioned there, it will be interesting to arrive at this result 
from first principles and check it in further models. However it is clear that no change in this 
contour prescription can correct the TBA ^-functions of || which will still have the wrong 
net overall change in g along the massive flows. 

We should also like to point out that the ^-functions we obtain to not satisfy the l g- 
theorem', but since this model is non- unitary that is perhaps no surprise. 

Another point to mention is that, in the massive case, Y(6,r) possesses square-root sin- 
gularities in the parameter r. These singularities lead to branch-points connecting the bulk 
vacuum energy to excited states [2(J, for which the Y-functions are described by so-called ex- 
cited state TBA equations [2C,22]. It would be nice to understand the role of these "excited" 



branches of Y in the framework of boundary theories, especially in the light of the result 6.1 
From the discussion in [^] it seems clear that in the massless case the excited branches Y n 
give the coefficients of the boundary state expanded in the eigenstates of the conformal field 
theory conserved quantities, i.e. if ( ip n | is the n-th excited state of the conformal field theory 
and | &(h) ) is the boundary state for the massless perturbation of the boundary <£, then (cf. 
eq. ( P3) ) 

(Vol <*>(/>)> = {<p\*(h)) =Y(f)(^\l) , (^„|*(fc)> =Y n (*£)(^|l) . 

Quite how this generalises to the massive model is not yet clear, but for the ground state we 
have given strong evidence for 



(nil 



Y(nr^) 



so it is reasonable to suppose that, by analytic continuation in A, the same should remain 
true in general, i.e. 

Some preliminary numerical and CPT checks support this conjecture. 

The TBA for ^-functions was used in a recent paper by Lesage et al. [|23[|, as part of an 
analysis of boundary flows in minimal models. In the light of the doubts raised above, it is 
possible that some of their conclusions for cases where the bulk as well as the boundary is 
perturbed should now be re-examined. However, since we have deliberately focussed on just 
one model in this paper, we will leave these questions for future work. 

Al. Zamolodchikov has recently found the exact correspondence between the UV and IR 
parameters in the Sine-Gordon model based on calculations of the extensive term f a from the 
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TBA ^-functions and from exact results for the expectation value of the boundary field \24\. 
Since the Yang-Lee model can be viewed as one special case of the Sine-Gordon model, his 
results should include our result relating h, M and b, summarised in equations (|4.4| ) and 
( 4.33; ), as a special case. So far, his results are not directly comparable with ours, but they 
are of the same functional form, only h c not being identified. A complete comparison between 
his results and ours would require the relation between the normalisations of the fields in the 
two approaches, and this has yet to be found. 
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